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ABSTRACT 

The very massive star system 77 Carinae exhibits regular 5.54-year (2024-day) period disrup- 
tive events in wavebands ranging from the radio to X-ray. There is a growing consensus that 
these events likely stem from periastron passage of an (as yet) unseen companion in a highly 
eccentric (e ^ 0.9) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hy- 
drodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and 
applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer 
(RXTE). By providing a global 3-D model of the phase variation of the density of the interact- 
ing winds, the simulations allow computation of the associated variation in X-ray absorption, 
presumed here to originate from near the apex of the wind-wind interaction cone. We find that 
the observed RXTE light curve can be readily fit if the observer's line of sight is within this 
cone along the general direction of apastron. Specifically, the data are well fit by an assumed 
inclination i = 45° for the orbit's polar axis, which is thus consistent with orbital angular 
momentum being along the inferred polar axis of the Homunculus nebula. The fits also con- 
strain the position angle that an orbital-plane projection makes with the apastron side of the 
semi-major axis, strongly excluding positions i/i < 9° along or to the retrograde side of the 
axis, with the best fit position given hy (f) — 27° . Overall the results demonstrate the utility 
of a fully 3-D dynamical model for constraining the geometric and physical properties of this 
complex colliding-wind binary system. 
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1 INTRODUCTION 

77 Carinae is one of most remarkable star systems in the galaxy. Its 
extreme luminosity, estimated today at some 5 x 10® Lq, implies a 
massive (M > lOOM©) primary star very close to the Eddington 
limit. One of the most extreme examples of the class of Luminous 
Blue Variable (LBV) stars, its historical light curve shows irregular 
brightenings, the greatest of which occurred in the 1840's, when its 
luminosity is estimated to have approached 25 x 10® Lq. This was 
accompanied by the ejection of some 10-20 A^q, forming what is 
seen today as the bipolar Homunculus nebula. In general, 77 Cari- 
nae is a key object in our understanding of the formation and evolu- 
tion o f extremely massive stars (see, e.g.. [Davidson & Humphrevsl 
fl997h . 
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An important advance in the observational study of 77 Cari- 
nae came from the identification of periodic, near-IR varia- 
tions (W hitelock et al . 1994; Damineli 1996) that are stable over 
many decades , along with co rrelated variability i n the radio 
( iDuncan et al.|[l995.) and X-ray dCorcoran et al.|[l99^ wavebands. 
The variability is especially dramatic in the 2-10 keV X-ray band, 
where the spatially unresolved X-ray flux drops by about a fac- 
tor of 100 for 3 months, as shown by daily monitoring with the 
Rossi X-Ray Timing Explorer (RX TE) during the X-ray minimum 
of 1997-1998 jlshibashi et ai][T999i), an d again during the 2003 X- 
ray minimum ( Hamaguchi et al.l 2007h . The top panel of figure [T] 
compares the RXTE lightcurve vs. phase over the two full orbital 
cycles for 1996-2001 and 2002-2007 (Corcoran 2005). 

Analysis of this X-ray emission and light curve has pro- 
vided important clues about the likely general nature of the sys- 
tem. First, the relative hardness of the X-rays suggests they must 
originate from the post-shock regions of a relatively fast wind 
(~ 3000 km s~^) from an otherwise unseen companion star, con- 
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Figure 1. Top: RXTE light curve vs. orbital phase, comparing data from 
the first (black) and second (red) orbital periods. Bottom: Comparison of 
the initial orbit RXTE light-curve (black) with our best-fit model (red), with 
observer position angles i = 45° and </> = 27°, as marked by "X" in the 
X contour plots in figure |4] For both comparisons, the lower boxes plot the 
square residuals of the fit, normalized so that the sum equals 1 . 



fined by the much den ser, but slower 500 — 800kms~^) 
wind from the primary (" Pittard et aT]|l998l : ICorcoran et alj|200ll : 
[Pittard & Corcoran 2002). The sharpness of the ingress and egress 
suggests moreover that the X-ray emission source must be rela- 
tively compact, probably originating mostly just inside the stag- 
nation point of the wind-wind shock cone, along the line between 
the stars. And given the very high density of the primary wind, the 
detection of X-rays during most of the period suggests an observer 
perspective that looks through a relatively transpare nt cavity carve d 
out by the relatively low-density secondary wind ( Corcoranll2005l) . 

A key hindrance to moving beyond this general picture has 
been the lack of a 3-D hydrodynamical wind-interaction model that 
fully accounts for the orbital motion, which can be especially im- 
portant for the sharp variations near periastron. The present paper 
applies Smoothed Particle Hydrodynamics (SPH) simulations to 
provide such a 3-D model throughout the full elliptical orbit of the 
binary components dOkazaki et al.l i2008). For simplicity, the ini- 
tial simulations here assume isothermal flow with a fixed common 
temperature for both winds. As such they do not directly model the 
shock-heated gas that is the cause of the X-ray emission. But the 
simulations do provide a fully 3-D, time-dependent description of 
the relatively cool material that is the source of X-ray absorption. 
By assuming a simple point-source model for the X-ray emission, 
located just within the head of the wind-wind shock interaction 
front (see figure the model allows computation of the phase- 
variable X-ray attenuation, and thus X-ray light curve, for any as- 
sumed observer position. As detailed below, with quite nominal bi- 
nary wind parameters adopted from previous analyses, the overall 
model, once adjusted to an optimal viewing angle, reproduces the 
observed RXTE light curve remarkably well (see figure[TJ. 



Table 1. Stellar, wind, and orbital parameters 



Parameters 


■q Car A 


ri CarB 


Mass (Mq) 


90 


30 


Radius (Rq) 


90 


30 


Mass loss rate (A/q yr~^) 


2.5 X 10"'' 


10-5 


Wind velocity (km s ^ ^ ) 


500 


3,000 


Wind temperature (K) 


3.5 X lO" 


3.5 X 10* 


Orbital period P (d) 


2,024 




Orbital eccentricity e 


0.9 




Semi-major axis a (AU) 


15.4 





2 MODEL SPECIFICATIONS 

The simulations presented here were performe d with a 3-D SPH 
cod e based on a versi on originally developed bv lBenz et al ]l ll990h 
and lBateetaIHl995l) . Using a variable smoothing length, the SPH 
equations with the standard cubic-spline kernel are integrated with 
individual time steps for each particle. In the implementation here, 
the artificial viscosity parameters are aspH ~ 1 and /3sph = 2. 

The two winds are modeled by an ensemble of gas particles 
that are continuously ejected with a given outward velocity at a 
radius just outside each star, coasting from there without any net 
external forces, effectively assuming that gravity is canceled by ra- 
diative driving terms. Perhaps more significantly, the simulations 
also assume both winds to be isothermal, with a common "warm" 
temperature. (The specific temperature, set to be comparable to the 
stellar effective temperature T = 35, 000 K, has little effect on 
the flow dynamics or X-ray absorption.) This is a serious simpli- 
fication, made to bypass the need to resolve the complex cooling 
regions near the wind shocks, which is generally difficult in a 3-D 
model, particularly with an inherently viscous method like SPH. 
While this does allow a quite realistic account for the 3-D absorp- 
tion by radiatively cooled material, it means that the expected X-ray 
emission from shock heating must be added separately (see §4). 

In a standard xyz Cartesian coordinate system, we set the bi- 
nary orbit in the x-y plane, with the origin at the system centre 
of mass, and semi-major axis along the x-axis (see figure |2l(. The 
outer simulation boundary is set at a radial distance r — 10.5a 
from the origin, where a is the semi-major axis of the binary or- 
bit. Particles crossing this boundary are removed from the simu- 
lation. By convention, we define t = (and zero phase) to be at 
periastron passage. Table [T] summarizes the stellar, wind, and or- 
bital parameters, larg ely a dopted from those derived previously by 
ICorcoran et"d] j200lh and lffillier et al.l ( l200ll) . A key parameter for 
the global form of the wind interaction is the ratio r/ of wind mo- 
mentum Mv between the primary to secondary wind, which here 
has a value r) ~ 4.2. Simple ram pressure balance then implies 
that, for a binary separation D, the interface should be located at 
a distance d = + •^/rj) ~ 0.33 D from the secondary star 

jStevens et alJll992l ; ICantoetai]| 19961) . 



3 PHASE VARIATION OF COLLIDING WINDS 

Figure [2] illustrates SPH simulation results for the density at 4 
phases from near periastron (left) to near apastron (right), plotted 
in both the orbital plane (x-y; top row), and the perpendicular plane 
through the orbital and major axes (x-z\ bottom row). Although in- 
stabilities in the wind-wind interaction lead to substantial stochas- 
tic variations and clumping, one can still see quite vividly how the 
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Figure 2. Snapshots of 3-D SPH simulation of wind-wind collision at times that, from left to right, are -1, +48, +202, and +1,011 days from periastron. 
The color scale shows the density (on a logarithmic scale with cgs units) in the x-y orbital plane (top) and in the x-z perpendicular plane containing the 
orbital and major axes (bottom). The main figures are for a square region ±10a about the system centre of mass, while the lower-left insets show a factor 
two magnification of the inner ±2a; this shows more clearly the interaction front where the X-ray source is assumed located, marked by an X between the 
brighter and dimmer spots that represent the primary and secondary stars, within the apex of the lower density wind. Annotations give the time (in days) from 
periastron passage and the number of particles (A'l, N2) in the primary and secondary winds. 

lower-density, faster wind from the secondary carves out a cavity in 
the higher-density, slower wind from the primary. Throughout most 
of the period centered around apastron, this cavity has a relatively 
simple, 2-D axisymmetric, conical form similar to the apastron 
snapshot at f = 1011 d, with a fixed opening ha lf-angle a ~ 60° . 
This is in good agreem ent with 2-D analytic ( Canto et"al]|l996l ) 
and numerical models dPittard & Corc oran 200^ that ignore or- 
bital motion, which near apastron is indeed small (ca. 20 kms~^) 
compared to the flow speed of either wind. 

But near periastron, the faster variation, closer separation, and 
higher orbital speed (up to ~ 360 kms^^) all work to distort the 
structure. In the approach up to periastron, the 2-D interface first 
starts to bend, but then, as the secondary whips around the oppo- 
site side of the primary, the secondary wind cavity becomes fully 
enshrouded by the denser, primary wind. Over time, the segment of 
this shell expanding toward apastron dissipates, and the nearly 2-D 
axisymmetric structure is again recovered. 





Figure 3. Schematic to illustrate observer position defined by inclination 
angle i to the orbital axis z, and by equatorial projection angle (/> relative to 
major axis x. The background orthogonal planes represent slices of density 
from our 3-D SPH simulation, shown here near apastron {t = 1011 d) for 
both the x-y orbital plane and the y-z plane orthogal to the major axis. The 
X-ray source is assumed to be just inside the head of the bow front, at a 
distance 0.25D from the secondary, where D is the current binary separa- 
tion. The observer direction shown is roughly consistent with that inferred 
by the best-fit model to the RXTE lightcurve. Note that, in the conventional 
notation of binary orbits, using a reference plane perpendicular to the line 
of sight, the "argument of periapsis" uj = 270° — </>. 



4 MODELING THE RXTE LIGHT CURVE FOR CAR 

To illustrate the diagnostic potential of this 3-D SPH simulation, 
we now use it to model the X-ray light curve observed by RXTE. 
The solid black curve in figure [T] shows this light curve for the 
years 1996-2007, covering the t wo initial full pe riods spanning both 
the 1998 and 2003.5 minima ( ICorcoranll200^ . While the sharp- 
ness of the drop to these minima seems suggestive of an eclipse- 
like event, the overall asymmetry does not fit the normal form 
of a stellar eclipse. The pre-event rise can likely be attributed to 
the 1/D scaling of the shock emission with the declining binary 
separation distance D toward periastron. But it is been a subject 
of debate whether the sharp drop and lack of a symmetric post- 
event peak reflects some kind of quenching of the X-ray emission 
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Figure 4. Le/f; Contours of reduced chi-squared, /{N — 2), for deviations between RXTE lightcurve for orbit 1 and models with observer inclination 
angle i and prograde angle <f>, plotted for a model grid with 9° increments in both i and 4>. The large "X" denotes the best-fit model plotted in figure[T] Right: 
Mosaic of lightcurves for inclinations and azimuths denoted by dots within the box bracketing the broad valley on left. The curves compare model (red) vs. 
orbit 1 RXTE data (black), with the difference shaded in cyan. The box in each panel gives the reduced chi-squared value. The red box highlights the best-fit 
model, withi = 45°, </> = 27°, and x^/{A^ - 2) = 0.51. 



jHamaguchi et alj2007h . or is mainly just due to variations in X-ray 
absorption. 

To explore the latter possibility, we combine the variable ab- 
sorption column derived from the SPH simulations with a simple 
point-source model for the X-ray emission. The strongest X-ray 
emission is expected to come from the shock of the faster secondary 
wind in the region just within the head of the wind- wind interaction 
front. In terms of the binary separation D, that interaction front is 
itself a distance d — D/{1 + ^) — 0.33 D from the secondary. 
Our model thus assumes a point source of X-ray emission located 
along the line of separation at a fixed fractional distance from the 
secondary, given by dx = fx D, where fx < 0.33. We have ex- 
plored models with fx — 0.20, 0.25, and 0.30, but since the results 
are all qualitatively similar, we focus here just on the intermediate 
case with fx = 0.25. 

Following the expected scaling for emission by adi- 
abatic shocks in wind-wind collisions tetevens et all 1 19921 : 
IPittard & Corcoranll2002) . we assume the phase variation of the X- 
ray source brightness varies with the inverse of the current stellar 
separation, Lx ~ 1/D. Defining then the time- variable mass col- 
umn depth from the X-ray source to the observer as rrioit), the 
model X-ray light curve takes the form. 



nod(i) = 



A 



+ B, 



(1) 



where A and B are normalization constants fixed to match the ob- 
served X-ray counts respectively at apastron and post-periastron 
minimum. Assuming a characteristic bound-free opacity k ~ 
5cm^ (see, e.g., fig. 5 of Antokhin et al. 2004) for the rele- 
vant RXTE energy band (2-10 keV), we then compute the phase 
variation of absorption from this X-ray source to trial observers 
over wide range of position angles. As illustrated in figure |2] this 
observer position is defined by the inclination i to the orbital axis, 



and by an orbital plane projection that makes a prograde direction 
angle with the +x axis direction toward apastron. 



5 VARYING OBSERVER POSITION FOR BEST FIT 

We have computed a grid of model X-ray light curves Lx^modit) 
vs. temporal phase t for a full range of observer's position angles i 
and (j), varying both in increments of 9° . The lower panel of figure[T] 
compares the first-orbit RXTE light curve with the resulting best- 
fit model, for which i — 45° and = 27°. The agreement is as 
good or somewhat better than the internal agreement between the 
first and second orbit cycles of RXTE observations, as shown by 
the black vs. red curves in the upper panel. In fact, the random vari- 
ations in model X-rays during the general rise before periastron ap- 
pear to be statistically quite similar to RXTE variations during this 
phase, though of course the random nature means they don't match 
in detail. In the model, these variations arise entirely from changes 
in absorption due to clumping in the wind interaction region of the 
SPH simulation, suggesting then that the observed variation might 
likewise be due to clump absorptions rather than, e.g., increases in 
the tempe rature or emission meas ure of the shock X-ray emitting 
region (cf. lHamaguchietal.ll2007h . 

For each observer position we quantify the level of agreement 
with the first RXTE orbit by the usual statistical measure of merit. 



X 



-y 



Lx ,mod 



1 



(2) 



where we have assumed the data can all be characterized by a com- 
mon /racftona/ mean-square deviation cr^. For the RXTE data, the 
contribution from measurment error is relatively unimportant com- 
pared with the inherent, apparently random variations in the ob- 
served X-rays, e.g. perhaps due to wind clumping. Moreover, in 
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figure[T]the comparison of the RXTE light curves for successive or- 
bital periods shows a systematic change, indicating a cycle-to-cycle 
variation that is not accounted for in our basic model. We thus esti- 
mate the inherent deviation by computing the averge mean-square 
deviation between each of these first two observation cycles, 

2 _1_ Lx^obsjt — P)i ^ 

where Lx.obs {t — P)i represents data from the second orbit shifted 
back by one period P = 2024 d and interpolated onto the data 
times of the first cycle. Application of this procedure for the RXTE 
data yields an estimated relative rms error, a = 0.20. 

Figure |4] plots contours of the reduced chi-squared, Xred = 
/{N — 2), for the most relevant subset of our model grid, with 
azimuth spanning a 90° range from just retrograde to strongly pro- 
grade of the major axis (—9° < (j) < 81°). Noting the overall 
north-south symmetry, the inclination spans the full range of just 
the northern hemisphere, < (j> < 90° . The formal best-fit model, 
marked with an "X", has observer position angles i — 45° and 
4> = +27°, with a Xred ~ 0.51 that is quite significantly below 
the unit value required for a good fit. (This suggests our derived 
a — 0.2 may be about a factor %/2 overestimate.) 

The contours also help identify the allowed range of view- 
ing angles, though this can be diffic ult to quantify rigorously. A 
common approach dPress et al is to define the difference in 

chi-squared relative to the best-fit model, which here gives 

Ax' = x\h 0) - xlmn = 479 X {xlei - 0.51) , (4) 

where 479 — N ~ 2 represents the number (A*' = 481) of data 
points in the first RXTE orbit, minus the two degrees of freedom 
(i, (f>) in the data fit. It turns out even models neighboring the best-fit 
have Ax' > 10, sometimes several tens or even in the hundreds; 
by formal statistics they would all be excluded at well above the 
99% confidence level. Taken at face value, this implies that, around 
the best-fit values i = 45° and (p — 27°, the range in both allowed 
viewing angles is less than the ±9° of the model grid. But this ap- 
proach very likely greatly overstates the real exclusion probability, 
given that we are fixing many model parameters about the orbit, 
winds, location of the X-ray source, etc. 

Nonetheless, even the reduced chi-squared contours do seem 
to strongly exclude azimuths with < 9° that are near or retro- 
grade of the major axis; likewise, inclinations near the orbital axis, 
i.e. with i < 36° , seem also excluded. On the other hand, viewing 
angles over the broad plateau within Xred ~ 1 (representing a dou- 
bling of the minimum) might still be allowed. The range in azimuth 
(9° < < 45°) and inchnation (36° < i < 90°) essentially just 
places the observer on the prograde side within the wind interaction 
cone of half-angle ~ 60° about the apastron side of the major axis 
X (see figure[2l(. 

The right panel of figure |4] compares lightcurves for viewing 
angles that bracket this allowed region. Comparison of the relative 
area of cyan shading between the observed and model curves sup- 
ports the view that the full range of inclination i = 45° — 81° 
around (p = 27° ± 9° give an acceptably good fit, while models 
outside this range do not. 

6 CONCLUSIONS AND FUTURE WORK 

The relative ease and natural way that the observed RXTE light 
curve is fit by this 3-D absorption plus point-source-emission 



model provides good evidence for the basic validity of the over- 
all paradigm, the key features of which are: 

(i) A highly elliptical orbit with the observer viewing from the 
general direction of apastron and prograde of the semi-major axis, 
through a cavity carved out in the slower, denser, primary wind by 
the faster, less-dense, secondary wind. 

(ii) A relatively localized X-ray source located on the secondary 
side of the interaction front between the stars. 

(iii) The X-ray mininum arising from a "wind eclipse" of this 
localized source as the primary wind engulfs the secondary wind 
just after periastron. 

Note that the last point implies that "quenching" of the X- 
ray emission is not likely to be a dominant effect in causing the 
broad-band X-ray mini mum. On the other hand, recent analyses 
jHamaguchi et al.l2007h suggest that some sort of spectral variation 
of emission may be necessary to explain observed changes in the 
X-ray hardness. In future work, we plan to extend our analyses to 
include more realistic models of the energy dependence of both the 
emission and absorption, with a particular focus on explaining such 
spectral energy and hardness variations. 
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